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A general approach to viable modified f{R) gravity is developed in both the Jordan and the 
Einstein frames. A class of exponential, realistic modified gravities is introduced and investigated 
with care. Special focus is made on step-class models, most promising from the phenomenological 
viewpoint and which provide a natural way to classify all viable modified gravities. One- and 
two-steps models are explicitly considered, but the analysis is extensible to Ai'-step models. Both 
inflation in the early universe and the onset of recent accelerated expansion arise in these models in 
a natural, unified way. Moreover, it is demonstrated that models in this category easily pass all local 
tests, including stability of spherical body solution, non-violation of Newton's law, and generation 
of a very heavy positive mass for the additional scalar degree of freedom. 
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I. INTRODUCTION 



Modified gravity models constitute an interesting dynamical alternative to the ACDM cosmology in that they are 
able to describe with success the current acceleration in the expansion of our Universe, the so called dark energy 
epoch. Moreover, modified F{R) gravity (for a review, see e.g. [l|) has undergone many studies which conclude 
that this gravitational alternative to dark energy [3, is able to pass the solar system tests. The investigation 
of cosmic acceleration as well as the study of the cosmological properties of F{R) models has been carried out in 
Refs. flSiiiia. 

Recently the im por tance of those models was reassessed, namely with the appearance of the so-called 'viable' F{R) 
models [i, [IS [nl, [12. Those are models which do satisfy the cosmological as well as the local gravity constraints, 
which had caused a number of problems to some of the first-generation theories of that kind. The final aim of all these 
phenomenological models is to describe a segment as large as possible of the whole history of our universe, as well as 
to recover all local predictions of Einstein's gravity, which have been verified experimentally to very good accuracy, 
at the solar system scale. 

Let us recall that, in general (see e.g. for a review), the total action for the modified gravitational models reads 

S^^J d^x^ [R + /(i?)] + . (1) 

Here /(i?) is a suitable function, which defines the modified gravitational part of the model. The general equation of 
motion in F{R) = R + f{R) gravity with matter is given by 

^g^uFiR) - i?^.i^'(i?) - g^.DF'iR) + V^V^F'iR) = -yT(,„)^, , (2) 

where T(^m)ni' is the matter energy-momentum tensor. 

In this paper we investigate two classes of 'viable' modified gravitational models what means, roughly speaking, 
they have to incorporate the vanishing (or fast decrease) of the cosmological constant in the flat (i? 0) limit, and 
must exhibit a suitable constant asymptotic behavior for large values of R. A huge family of these models, which 
we will term first class — and to whom almost all of the models proposed in the literature belong — can be viewed 
as containing all possible smooth versions of the following sharp step-function model. To discuss this toy model, at 
the distribution level, will prove to be very useful in order to grasp the essential features that all models in this large 
family are bound to satisfy. In other words, to extract the general properties of the whole family in a rather simple 
fashion (which will involve, of course, precise distribution calculus). 

This simple model reads 
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FIG. 1: Typical behavior of f{R) in the one-step model). 

where 9{R — Rq) is Heaviside's step distribution. Models in this class are characterized by the existence of one or 
more transition scalar curvatures, an example being Rq in the above toy model (but there can be more, as we will 
later see). 

The other class of modified gravitational models that has been considered contains a sort of 'switching on' of the 
cosmological constant as a function of the scalar curvature R. A simplest version of this kind reads 

fiR) - 2Aeff(e-''« - 1) . (4) 

Here the transition is smooth. The two above models may be combined in a natural way, if one is also interested 
in the phenomenological description of the inflationary epoch. For example, a two-steps model may be the smooth 
version of 

f{R)^-2Aoe{R-Ro) -2Aje{R~Ri), (5) 

with i?o << Ri, the latter being the inflation scale curvature. 

The typical, smooth behavior of f{R) associated with the one- and two-step models is given, in the smooth case, 
in Figs. [1] and [2l respectively. The main problem associated with these sharp models is the appearance of possible 
antigravity regime in a region around the transition point and antigravity in a past epoch, what is not phenomeno- 
logically acceptable. On the other hand, an analytical study of these models can be easily carried out, as discussed 
in the Appendix. 

The existence of viable (or "chameleon" ) f{R) theories with a phase of early-time inflation is not known to us from 
the literature. The fact that we are able to provide several classes of models of this kind that are consistent also with 
the late-time accelerated expansion is thus a novelty, worth to be remarked. 

The paper is organized as follows. Next section reviews f{R) gravity in the physical, Jordan frame. Equations of 
motion are presented and the solar system tests (absence of a tachyon, stability of the spherical body solution, and 
non- violation of Newton law) are discussed. In Sect. Ill we give a number of viable modified gravities which may lead 
to the unification of early-time inflation with late-time acceleration and satisfy the solar system tests. Some properties 
of such viable modified gravity are discussed in detail. Sect. IV is devoted to the presentation of viable models via 
conformal transformation, as a kind of scalar-tensor theory which is mathematically equivalent to the original theory. 
In some examples, the explicit form of the scalar potential is derived. Corrections to Newton's law are also obtained 
for some of the realistic theories here considered. It turns out that these models do pass the stringent Newton law 
bounds, since the corrections to Newton's law in these cases turn out to be negligible. Some summary and outlook is 
given in the Discussions section. Finally, in an Appendix we show how to evaluate the positions of the corresponding 
de Sitter critical points, by considering the sharp version of the one- and two-step models, expressed in terms of the 
Heaviside and Dirac distributions. 

II. MODIFIED GRAVITY IN THE JORDAN FRAME 

Regarding the precise determination of the modifying term, f{R), we here revisit this issue in the Jordan frame 
(instead of the Einstein one). Let us recall the two sufficient conditions which often lead to realistic models (see, for 
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FIG. 2: Typical behavior of f{R) in the two-step model. 



example [§|) 

/(0) = 0, lim/(i?) = -a, (6) 

H — >Hi 

where a is a suitable curvature scale which represents an effective cosmological constant, being i?i >> Rq, with 
Rq > 0, the transition point. The condition /(O) = ensures the disappearance of the cosmological constant in the 
limit of flat space-time. 

By using these conditions, some models in this class are seen to be able to pass the local tests (with some extra 
bounds on the theory parameters) and are also capable to explain the observed recent acceleration of the universe 
expansion, provided that a = Aq = 2Hq, Hq being the Hubble constant at the epoch of reference. However, they do 
not incorporate early-time inflation, which comes into play at higher value of R. Thus, one might also reasonably 
require that [ll| 

/(0) = 0, \im f{R) = -{a + aj), (7) 

H — >H2 

where a/ >> a is associated with the inflation cosmological constant. A/, and where R2 >> Ri >> Rq, Ri being 
the corresponding transition large scalar curvature. 

Further restrictions, like small corrections to Newton's law and the stability of planet-like gravitational solutions 
need to be fulfilled too [|. All those can be also formulated in the mathematically equivalent Einstein frame. Here we 
present a short review of them in the Jordan frame which we consider as the physical one (see, Ref. @ for a discussion 
of the physical (non-) equivalence of the Einstein and Jordan frames). 

The starting point is the trace of the equations of motion, which is trivial in the Einstein theory but gives precious 
dynamical information in the modified gravitational models. It reads 

3V^f'{R) = R + 2f{R)-Rf'{R)~K^T. (8) 

The above trace equation can be interpreted as an equation of motion for the non trivial 'scalaron' f'{R) (since it is 
indeed associated with the corresponding scalar field in the other frame). For solutions with constant scalar curvature 
i?* , the scalaron field is constant and one obtains the following vacuum solution: 

i?»+2/(i?,)-i?,/'(i?,) = 0. (9) 

Furthermore, according to we can describe the degree of freedom associated with the scalaron by means of a 
scalar field x> defined by F'{R) = 1-1- f'{R) = e~^. If we consider a perturbation around the vacuum solution of 
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constant curvature i?, , given by i? = i?* + 5R, where 



then the equation of motion for the scalaron field is 

™^ " 3 l^^^ - j = - 6(l + /-(i?.) ^ ■ ^''^ 

As a result, in connection with the local and with the planetary tests, the following effective mass plays a very crucial 
role: 

If < 0, a tachyon appears and this leads to an instability. Even if NP > 0, when AP is small, it is SR 7^ at 
long ranges, which generates a large correction to Newton's law. As a result, has to be positive and very large in 
order to pass both the local and the astronomical tests. This stability condition can be also derived within QFT in 
de Sitter space-time (see, for instance, [ijl) . 

Concerning the matter instability (sl. llSl. [l6j. this might occur when the curvature is rather large, as on a planet, 

as compared with the average curvature of the universe R ~ (10~'^^ eV) . In order to arrive to a stability condition, 
we can start by noting that the scalaron equation can be rewritten in the form [l5| 



r(i?) (i+/mR 2{R+f{R)) _ . 
°^ + ' ~W{R) WW~ WW) • ^ ^ 



If we now consider a perturbation, 5R, of the Einstein gravity solution R = R^ = > 0, we obtain 

0^ {~d^ + U{Re))SR + C , (14) 

with the effective potential 



F""(i?e) F'"{R,)^\^ „ „ , i?e 



F'{Re)F"'iR,)Re F'{R,) , 2F{R,)F"' {R,) F"'{Re)Re 



3F"{Re)^ 3F"{Re) ?,F"{Ref ?,F"{ReY 



(15) 



If U{Re) is positive, then the perturbation 5R becomes exponentially large and the whole system becomes unstable. 
Thus, the matter stability condition is, in this case, 

[/(i?e)<0. (16) 

Coming back to the vacuum condition ([9]), we recall that, within the cosmological framework, it may be rederived 
making use of the dynamical system approach. This consists in rewriting the generalized Friedmann equations of 
the modified gravitational model in terms of a first-order differential system and looking for its critical points. For a 
modified gravity model, the associated dynamical system can be written as [71, 181] (a(t) being the expansion factor in 
a FRW flat space-time): 

a in a 

■J—np - [2(2 ~ nR)-3{w + l) + I -Qf-^p]^p, (17) 
a In a 

where w = ^ is the usual barotropic constant, being 
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with 



There exists yet another quantity 



which satisfies the constraint 



1 + m 



fjR) 

H{l + f'{R)) 



^, = - rJ::L.. ^ (20) 



np + Qp + ^p^i- (21) 

The critical points are solutions of the (algebraic) system 

= 2nFi'2 ~ ^Ir) + i^F - nR){i - ftp - ^p) , 

= [2{2- Qr) ^3{w + I) + I -Qf -^p]^p- (22) 

As an example of those, the de Sitter critical points are the ones in the invariant vacuum submanifold Qp = 0. 
These solutions read 

nR^2, R, = l2Hl , (23) 

and 

nF = l, R, = RJ'iR,.)-2f{R,). (24) 

The last equation coincides with Eq. It is a transcendental equation, for all our models, and can be solved only 
by iteration (see Appendix) or either numerically by other methods. 

The stability condition associated with the de Sitter critical point in this dynamical system framework can be 
investigated too, and reads @ 

'<««-) = 1^7W- '''' 



It coincides with the requirement that the effective mass ([T2|) be positive. In the matter-radiation sector, where fi^ 
is non- vanishing, other critical points may also exist. In the next section, we give explicit examples of exponential, 
viable modified gravity. 

III. EXAMPLES OF REALISTIC EXPONENTIAL MODIFIED GRAVITY 

After the general discussion above, we will here present some new viable f{R) models. We start with a most simple 
one 

fiR) = aier"'' - 1) . (26) 

Since /(O) = and f{R) —a for large R, conditions ^ are satisfied. Moreover, 

fiR) = -5ae~''^ , /"(i?) = b^ae"^'^ . (27) 

We have seen that in the discussion of the viability of modified gravitational models, the existence of vacuum constant 
curvature solutions plays a very crucial role, namely the existence of solutions of Eq. ([9]). With regard to the trivial 
fixed point i?* = 0, this model has the properties 

l + /'(0) = l-a6, f"{0) = ab\ (28) 

Thus, the effective mass for i?* = is 
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and then Minkowski space time is stable as soon as a6 < 1. Such condition is equivalent to 1 + /'(O) > 0. 

In order to investigate the existence of other fixed points, we first have to find the existence conditions and then 
make use of Newton's method or some of its variants (see the Appendix). It is easy to see that for the model (|26p . 
one has the critical point only if ab > 1, namely K'{0) > 1, where the function K{R) defined in Appendix. Then, 
one can construct an approximation procedure to solve Eq. ([9]), in terms of an iteration process, namely 

p _ p _ Rn - Rnf'jRn) + 2f{Rn) , , 

" l + /'(i?„)-i?„/"(i?„) ■ ^^"^ 

For a starting point, i? = large enough, f{R) is approximately constant and a few iterations give 

i?*,i~2a. (31) 

This critical point is stable and corresponds to the current acceleration in the universe expansion. This follows from 
that fact that the effective mass is 

- ^e^"" , (32) 

namely it is positive and large. However, since ab > 1, for very small i?, one has antigravity effects, namely 
1 + f'{R) < 0, as we shall see in Section IV. 

Matter instability can also be investigated. Eq. (|15p gives 

C/(i?,)^-^e2''«=, (33) 

which is negative, thus the matter stability condition is fulfilled. 

The model we have discussed so far does not exhibit a sharp transition curvature. But there are many models where 
one or more transitions of this kind appear. The Hu-Sawicki (HS) model 0] belongs to this one-step class family of 
models. A simple choice for a one-step model is in our case 



f{R) 

For R very small, we have 



1 + 6-''^° 



bR _ I 



m^~Y^^^ + oiR'), (35) 



while for suitable values of 6, one has the same behavior as in the HS model, where the continuous parameter b plays 
the role of the integer n in the above mentioned model. Higher values of b give rise to a sharper transition, occurring 
at Rq from very small values of f{R) towards a constant value —a. This model has an effective mass, evaluated at 
i?* = 0, which turns out to be negative, thus Minkowski space-time is unstable in this case, and it might happen that 
1 + f'{R) < around the transition. 

A simple modification of the above model which incorporates the inflationary era, namely the requirement ([7]), is 
a combination of the two models discussed so far, that is 



or, as a two-step model. 



bR _ 1 

fiR)=a{e-''^-l)-ar (36) 



gbi? _ I ^bR _ I 

fiR) ^ -a ^^^ ^ ^^^^ - ai ^ ^^^^ . (37) 



Again, /(O) ~ and, at the value R — Ri, there is a transition to a higher constant value —{a + aj) which can be 
related to inflation. 

We should note that f{R) in ([37| is a monotonically decreasing function. Then, when R is very large, f{R) tends 
to a constant value, which could correspond to the effective cosmological constant generating inflation. In order to 
describe the recent accelerating expansion of the universe, f{R) should remain almost constant, that is, f'{R) ~ 0, 
for sufficiently small values of R corresponding to the curvature of the present universe. 
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The main problem with this models is the appearance of antigravity in connection with the transition point when 
the function is sharp (see the discussion in the Appendix). The appearance of antigravity in the past, namely around 
i?/, is not acceptable, as we have already commented. 

A possible modification of the previous model is the following: 

fiR) = - 1) + cR^ ^^^_^^,^^ , (38) 

with N > 2 and c > 0. In this variant, similarly to the theory [l^, during the inflationary era at R > Rj , f{R), the 
model acquires also a power dependence on the scalar curvature, which may help to exit from the inflationary stage. 

As discussed in detail in the Appendix, for the sharp, theta models, besides the problem of antigravity, for Rq << a 
and Rj << aj, they posses, generically, two De Sitter critical points, one around the transition point i?* ~ and 
the other being 

-R*,2 ^ 2a . (39) 
We can also investigate the matter instability. For the two-step model (|37p . we now assume 

i?o < i? i?e < i?/ ■ (40) 

Then f{R) in p7)) can be approximated as 

f{R) ^-a{-l + {l + e-^"") e-^(«-«c)} - . (41) 

We may assume 

since bRj could be very large (see the argument around ([75)1 about antigravity). Then we find 

g6(«e--R.0) 

- ' 3a6Ml + e^^^o) ' ^''^ 

which is negative and there is no instability. 

We conclude this Section with a variant of the above model which facilitates the analytic computation and the 
discussion concerning antigravity. In fact, as a smoothed one-step function, we may consider 

f , fb{R-Ro)\ , fbRoW feKR-Ro)_i e''^" - 1\ 

f{R) = -a (^tanh ^ ' j + tanh ( ^ ) ) ^ ( ^^^^^T?:^ + 37?^ ) (44) 

When i? ^ 0, we find that 



f(R) ^ , . (45) 

' 2cosh2(^) ^ ' 



and thus /(O) — 0, as required. On the other hand, when R 



fiR) ^ -2Aeff = -a ( 1 + tanh ( ^ ) ) . (46) 



2 

li R^ Rq in the present universe, Aoff plays the role of the effective cosmological constant. We also obtain 

fiR) 7^ r , (47) 

2cosh2(^(^) ^ ' 

which has a minimum when R = Rq: 

nRo) = -Y- (48) 
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Then in order to avoid antigravity, we find 



0< l + /'(i?o) < 1 



ab 



(49) 



The model given by Eq. (|44| is able to describe late acceleration. In order to show that the de Sitter critical points 
exist, we can compute the function K{R) — Rf'{R) — 2f{R) of the Appendix and we have K'{R) — Rf"{R) — f'{R), 
K"{R) = Rf"'{R), 1 - ^'(0) = 1 + /'(O), where 



K'{R) = - 



bR tanh 



b{R-Ro) 



fiR), 



Thus, enforcing the absence of antigravity, one has 



1 + /'(O) > . 



(50) 



(51) 



namely, K'{0) < 1. In this case however K'{R) > for R > Rq, and it is negative for R < Rq. Therefore, we cannot 
use the argument of the Appendix. For the one step-model, however, one can actually live with antigravity in the 
future, thus in the sharp version, the analysis in the Appendix leads again in fact to the existence of two dS critical 
points. 

As a model which is able to describe both the inflation and the late acceleration epochs, we can consider the 
following two-step model: 



/(i?) = — ao ( tanh 
We now assume 

and 



bo (R - Ro) 



tanh 



boRo 



aj tanh 



bi (R - Ri) 



tanh 



hiRi 



Ri '> Ro , a/ » ao , bi <€. bo 



biRi > 1 



When i? or i? < i?o, Ri, f{R) behaves as 



fiR) 



aoba 



a lb I 



\^2cosh2(^) 2cosh2(^!^; 
and find /(O) = again. When R Rj, we find 



R . 



f{R) 



-2Ar 



ao 1 + tanh 



boRo 



— a/ 1 + tanh 



On the other hand, when Ro R Rj, we find 



1 + tanh 



boRo 



fiR) ^ -ao 
Here we have assumed (|54|) . We also find 

fiR) = - 



aibjR 



f ^iRi 



-2Ao = -ao 



—a/ 1 + tanh 



fbiRi 



\ 2 



1 + tanh 



boRa 



aobo 



otibi 



2coshMMZ^^ 



2 cosh^ 



which has two valleys when R ^ Ro or R ^ Rj . When R ^ Rq, we obtain 

a lb I 



f'{Ro) = -aobo - 



2 cosh^ 



2 / bi{R„-R.i} 



> —aibj — aobo 



(52) 

(53) 
(54) 

(55) 

(56) 
(57) 
(58) 



(59) 
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On the other hand, when R = Rj, we get 



f'(Ri) — —ajbi ^ > —ajbi — aobo ■ (60) 



Then, in order to avoid the antigravity period, we find 

aibi + aobo < 2 . (61) 

The existence of the de Sitter critical points in this two-step model is much more difficult to investigate. However, in 
order to get the acceleration of the Universe expansion it is sufficient that ujef / < — | . 

We now investigate the correction to the Newton's law and the matter instability issue. In the solar system domain, 
on or inside the earth, where R ^ Rq, f{R) in (|44| can be approximated by 

f{R) ~ -2Aeff + 2ae-''(-^--^") . (62) 

On the other hand, since i?o <C -R ^ Ri, by assuming Eq. ([5^ . f{R) in ([5^ could be also approximated by 

f{R) - -2Ao + 2ae-'"'(-"^-"°) , (63) 



which has the same expression, after having identified Aq = Acff and bo = b. Then, we may check the case of ((62)) 
only. 

We find that the effective mass has the following form 

which could be very large again, as in the last section, and the correction to Newton's law can be made negligible. 
We also find that U{Rb) in (fT5|) has the form 

U{Re) = (^2A + e-^(«-«°^ , (65) 

which could be negative, what would suppress any instability. 

Thus, we have here presented several realistic exponential models which naturallyunify the inflation with the dark 
energy epochs (with a radiation/matter dominance phase between, as in Refs. In addition, the Newton law 

is respected and all spherical body solutions (Earth, Sun, etc) are stable. 



IV. VIABLE MODELS IN THE EINSTEIN FRAME 



As is well known from previous studies, it is often quite convenient to go from the Jordan (physical) frame to 
the mathematically-equivalent Einstein frame description, where f{R) models become scalar-tensor theories with a 
suitable potential. In particular, corrections to Newton's law and the matter instability can be also investigated in 
the Einstein frame directly, where the relevant degrees of freedom are a new tensor metric and a scalar field. More 
specifically, concerning the infiation issue, the Einstein frame can indeed be very useful. Following e.g. reference 
we can introduce the auxiliary field A and rewrite the action ([T]) as 

S=-^J d^x^g {(1 + nA)) {R~A)+A + f{A)} . (66) 



From the equation of motion with respect to A, if f"{A) ^ 0, it follows that A = R. By using the conformal 
transformation g^i, — s- e'^g^^, with cr = — In (1 -I- /'(A)), we obtain the Einstein frame action [3|: 

- 1 [ d^,^r-AR-U^\\^^d Ad A~^ + I^\ 



^ J d^x^g (r - ^gP-dpad,a - Via] 



(67) 



Via) = e^g (e"-) - e^^F {g (e"-)) = ^ ^ • (68) 
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Here g {e~"') is given by solving a = — In (1 + /'(A)) = \nF'{A), as A = g{e~''). After the scale transformation 
9^iv ~^ gfiv is done, there appears a coupling of the scalar field a with matter. For example, if matter is a scalar 
field <I>, with mass M, whose action is given by 

^4> = IJ d^x^ {-g^^'^d^^d,^ - M2$2) ^ (69) 

then there appears a coupling with a (in this Einstein frame): 

S^E^\j d^xV^ (-e'^5^''9^$c),$ - M2e2-$2) . (70) 

The strength of the coupling is of the same order as that of the gravitational coupling, k. Unless the mass corresponding 
to cr, which is defined by 



,2 



, _ 1 d^Vja) _ 1 J A AFjA) 1 , 

" ~ 2 da^ 2\f'{A) {F'iAjf F"{A) ' ' ^ ' 



is big, there will appear a large correction to the Newton law. Newton's law has been investigated in the solar system, 
as well as on earth, where the curvature is much larger than Rq. For the model (j26p . we find 

bR 

which is positive. As we will discuss soon f(fM l) - (|78p ). we find 1/6 <^ Rq <^ R, and therefore bR ^ 1, which tell us 
that could be very large and the correction to the Newton law would be very small. For the model ([57]) . we find 

^ 2ab^l + e-bRo) ' ^^^^ 

which could be very large again, and the correction to the Newton law correspondingly very small. 
Eq. (IMl) also tells that if 



1 + /'(A)<0, (74) 
then antigravity could appear, since the effective gravitational constant is given by 

'^off ^ TTJiA) ■ ^^^^ 

In order to avoid it, the condition (|75p must be satisfied, at least until present, all the way since the beginning of the 
universe. Some remark is in order. In the antigravity region, there is no evolution of the universe with a flat spatial 
part. In the usual Einstein gravity, we have the FRW equation, {3/k^)H^ = p, but in the antigravity region, since 
the sign of changes, as —k^, we get —{3/k^)H'^ — p. Since the Ihs of this equation is always negative and 

the rhs is always positive, there is no solution, what shows that there is no time-evolution of the universe in this case. 
We should also note that, even if 1 + f'{A) is negative, the conformal transformation itself can still be well defined, if 
we use the absolute value of 1 + f'{A), that is, |1 + f'{A)\g^„ g^i,. In that case, however, in the obtained Einstein 
frame action, the sign of the scalar curvature R becomes negative, that is, antigravity again appears. Moreover, the 
initial value problem which is formulated in f{R) gravity via conformal transformation 'l3'| is not well defined. This 
is the reason why we avoid the consideration of antigravity regimes. 
For the simple model the condition ([7^ reads 

1 - abe-''^ > . (76) 

Since the scalar curvature R in the past universe could be larger than the curvature Rq in the present universe, we 
find 1 - afoe"''^ > 1 - abe'''^" . Therefore, if 

1 - abe-''^° > (77) 
is satisfied, the condition ([71]) can be satisfied too. Eq. (|77p tells us that 

a6 < 1 or i < i?o . (78) 
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In order that f{R) can play the role of an effective cosmological constant for the present universe, the second condition 
1/b Rq should be preferred. The situation is not much changed in the two-step model p7p. and the antigravity 
condition tells us that 1/6 <C Rq, again. 

We can now investigate the regions which reproduce realistic models. For the simple model (|H), conditions ([77]) or 
([78)) . which avoid antigravity in the history of the universe, can be expressed as 

i < i?o (10-33 eV)^ . (79) 

On the other hand, the condition that the correction to Newton's law should be small is that nia, given by ([7^ . must 
be large. Since a in (|72p plays the role of the effective cosmological constant in the present universe, we have 

a - i?o (lO-^^eV)^ . (80) 

In the solar system, we find R - lO^^^ eV^ Even if we choose 1/6 ~ i?o - (lO^^^ ^ that ml - IQi^ooo eV^ 

which is ultimately heavy. Then, there will not be any appreciable correction to the Newton law. In the air on earth, 
R - 10-^°eV^ and even if we choose 1/6 ^ Rq ^ (lO-^^eV)^ again, we find that ml - ioio,ooo,ooo,ooo Then, 
the correction to Newton's law is never observed in such model. 

For the model (l37|) . since _Ro <C i? <C i?/ in the solar system or on the earth, f{R) can be approximated by ([4T|) . 
Then the effective gravitational constant could be given by 

1 _ 1 / a/6 \ 

The mass of the scalar field a is given by ([75]) . which is very large again, that is, ml ^ 10^'°°" eV^ in the solar system 
and ml ~ j^qio, 000,000,000 ^^le air surrounding the earth, and therefore the correction to the Newton law is 

negligibly small, either. 

There is a technical point which deserves more careful considerations. It is that A = R has to be expressed as a 
function of a by solving the equation 

fiA) =e---l (82) 

and this can be explicitly done for the simplest cases only. For example, in Ref. Q a class of models defined by means 
of the function 

1 + C2(/t/TO ) 

has been proposed. Here ci, C2 are arbitrary dimensionless constants, while m has the dimension of mass. This model 
yields an effective cosmological constant which generates the late-time accelerated expansion. For such class of models, 
Eq. ([5^ reduces to an algebraic equation of order 2n, which can be explicitly solved for n — 1 and n — 2. In the 
simplest case, n — 1, one easily gets 



A± = — 

C2 



e'^ - 1 



C2 > . (84) 



Via) = e^il-e^) A-e'-f{A) 



C2 



(85) 



where the positive solution has been chosen. For n — 2 the potential assumes a quite complicated form, which is 
practically useless. 

A simple modification of the model (|83p is the following ^iTl: 



2ci(i?/m2)" -hC3 



l + C2{R/m'^y 
which for n = 1 and C2 > gives rise to the potential 



n > 1. 



(86) 



Via) 



C2 



(ci -f l)e" - 2e"/V(ci-C2C3)(e--l)- 1 



(87) 
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Now, we go back to the models that we have considered above. For some of them we can give an expUcit form for 
the potential. We start with our first model In such case, Eq. ([5^ is a simple transcendental equation which 
gives rise to 



-I- 



ab ^ 



V(.) . - 



l + {ab- 1)6" + {e" ~ 1) In 



1 



ab 



(89) 



where cr > is understood. 

Also for the one-step model (|34|) . equation ((82)) becomes a transcendental one, but it can be solved and we can 
eventually write the potential in the form 



where A = A{a) is given by 



(90) 



A, 



log ■ 



27(6"^ - 1) 



i-^ii^±Ji-fci 

M V M 



H = a/3(l + 7) 



where 7 — e~^^° . 

Also for the two-step model p7l) . which includes inflation, one can obtain an exact expression for the potential 
V{a) but, since a/ 3> a, such expression reduces to the latter above, with the replacement a aj. The explicit 
expressions of the scalar potential in the equivalent, scalar-tensor theory can be actually very useful in the study of 
the PPN-regime of modified gravity, in that of the stellar evolution equations, and also in some related questions. 



V. DISCUSSION AND CONCLUSIONS 



In this paper, a general approach to viable modified gravity has been developed in both the Jordan and the Einstein 
frames. We have focussed on the so-called step-class models mainly, since they seem to be most promising from the 
phenomenological viewpoint and, at the same time, they provide a natural possibility to classify all viable modified 
gravities. We have explicitly presented the cases of one- and two-step models, but a similar analysis can be extended 
to the case of an A^-step model, with A^ being finite or countably infinite. No additional problems are expected to 
appear and the models can be adjusted, provided one can always find smooth solutions interpolating between the de 
Sitter solutions (what seems at this point a reasonable possibility) , to repeat at each stage the same kind of de Sitter 
transition. We can thus obtain multi-step models which may lead to multiple infiation and multiple acceleration, in 
a way clearly reminiscent of braneworld infiation. 

This looks quite promising, with the added bonus that the model's construction is rather simple, as we have here 
shown explicitly. All the time, as a guide for an accurate analysis, use has been made of the simple but efficient 
tools provided by the corresponding toy model constructed with sharp distributions, a new technique that we have 
here introduced too. It is to be remarked that, for the infinite-step models, one can naturally expect to construct the 
classical gravity analog of the stringy landscape realizations, as in the classical ideal fluid model [l7j . 

For the model ([5^ . both inflation in the early universe and the recent accelerated expansion could be understood 
in those models in a unified way. If we start with large curvature, f{R) becomes almost constant, as in (j56p . and 
plays the role of the effective cosmological constant, which would generate inflation. For a successful exit from the 
inflationary epoch we may need, in the end, more (say small non-local or small i?") terms. When curvature becomes 
smaller, matter could dominate, what would indeed lower the curvature values. Then, when the curvature R becomes 
small enough and Rq <^ R Rj , f{R) becomes again an almost constant function, and plays the role of the small 
cosmological constant which generates the accelerated expansion of the universe, that started in the recent past. 
Moreover, the model naturally passes all local tests and can be considered as a true viable alternative to General 
Relativity. Some remark is however in order. On general grounds, one is dealing here with a highly non-linear system 
and one should investigate all possible critical points thereof (including other time-dependent cosmologies), within 
the dynamical approach method. Of course, the existence of other critical points is possible; anyhow, for viable f{R) 
models, to find them is not a simple task, and we have here restricted our effort to the investigation of the dS critical 
points. With regard to the stability of these points, the one associated with inflation should be unstable. In this way. 
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the exit from inflation could be achieved in a quite natural way. In particular, for instance, this is in fact the case for 
the two step model with the term discussed in the Appendix. 

In conclusion, a class of exponential, realistic modified gravities have been here introduced and investigated with 
care. Some of these models ultimately lead to the unification of the inflationary epoch with the late-time accelerating 
epoch, under quite simple and rather natural conditions. What remains to be done is to study those models in further 
quantitative detail, by comparing their predictions with the accurate astrophysical data coming from ongoing and 
proposed sky observations. It is expected that this can be done rather soon, having in mind the possibility to slightly 
modify the early universe features of the theories here introduced, while still preserving all of their nice, realistic 
current universe properties, as we have shown above. 
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APPENDIX A 



In this Appendix, we will study how to evaluate — or at least get some information — on the number and positions 
of the corresponding de Sitter critical points. In order to grasp the general behavior of the larger family of models (as 
already discussed in Sects. I and II), let us start by considering the sharp (mathematically very clean albeit physically 
unrealistic) version of the one-step models, expressed in terms of Heaviside and Dirac distributions, namely 

,f(R)^-ae{R- Ro), f(R) = -aS{R-Ro), a>0. (Al) 

This simple, idealized model leads to antigravity, since 1 + f'{R) is obviously always negative at the transition. In 
practical terms, this means that the sharper is the smoothing of the step-function, the harder one will be involved in 
the antigravity problem. For the one-step model, antigravity could be arranged to happen in the future. In fact, the 
equation whose solutions are the de Sitter critical points becomes 

R= -a Ra S{R - Rq) + 2a9{R - Rq) . (A2) 

This is an equation in the distribution theory sense and requires an appropriate treatment. To start, there is the 
trivial solution R — 0. li Ry^ Rq, the only solution is R — 2a. However, eventually we have to deal with a non-ideal, 
physical situation and we must consider not the sharp but the smoothed version of the delta and theta distribution. 
For the delta, we may consider its support to be contained in the interval —e + Ro, Ro + £■ 

We can get information about the other fixed point by arguing as follows. Integrating the above equation from the 
value R = Rq — e to the value R{) + e: 

1 {ARos) = -aRo + laRoe . (A3) 
As a result, we obtain the consistency condition for e 

Typically, we have Rq « a. In this case, we obtain 

^<i?,i<^ i?.,.^2a. (A5) 

Thus, with a sharp f{R) function, one has two de Sitter solutions when Rq « a. antigravity effects can be confined 
around (namely, in the future) and the current acceleration is represented by the second solution i?*^2 — 2a. 

For a two-step sharp model, however, this solution to the problem is not acceptable, because we cannot allow for 
antigravity in past epochs, as discussed in Sect. IV. The way out of this is to consider, for example a sufficiently non- 
sharp smoothing of the theta functions, but in this case, the above analysis is not longer valid. Another possibility 
— which can still make use of sharp theta functions — is the following modification of the two-step model, by a power 
of the curvature: 

f{R) = -a0{R- Ro) + f3R^9{R- Ri), a,/?>0. (A6) 
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In this case, near the second transition point one does not have, by construction, antigravity effects. The above 
analysis can now be appUed again, with the result that, integrating from Rj, iV > 2, TV = 3, etc., one gets 

i?* ~ 2a + (3{N - 2)i?f . (A7) 

This is an algebraic equation of N-th order, whose solutions can be easily investigated, in any specific situation. For 
example, for TV = 3 and 2a negligible, one gets the approximate value 

. (A8) 

In contrast, in the non sharp case it is not so easy to find sufficient conditions for the existence of the de Sitter 
critical points. With regard to this issue, let us recall the following theorem, which may be useful in a direct, numerical 
computation of the critical points in a physically realistic setting. 

Theorem: Given a twice differentiable function, G{x), defined in the real interval [a,b] and such that G'{x) is non- 
vanishing in this interval, and 

\G"{x)G{x)\<\G'{xf\, (A9) 
then the zeroes of G{x) are obtained by the recursive formula (Newton's tangent method) 

Xn+l = Xn- 7^77 — T . (AlO) 

Gr (Xn) 

This theorem yields a contraction mapping on the complete metric, in our case in the interval [a, b]. In fact the zeroes 
of G{x) are the fixed points of the function K{x) — x ~ g^(x) ' ^^^^ function is a contraction mapping, as far 
as < 1. Since K'{x) — ^^-^^f^, one gets the stated result. The recursive relation Xn+i — K{xn) leads to 

Eq. (Km . 

In the cases we consider here, this result is, however, not easy to implement. Alternatively, one can proceed as 
follows. Let us write 

K{x) ^ xf'{x) -2f{x) , K{0)^0, K{x) -2f{x) , x»0. (All) 

The fixed points x = K{x) are the de Sitter critical points. Let us suppose that K'{x) > 0, for every a; > 0. Then it 
follows that: (i) if K'{0) > 1 then there exists a fixed point; (ii) if K'(0) < 1 then there are no fixed points. 
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